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1. Introduction 

There is an extensive experimental effort aiming at a precise determination of neutrino oscil- 
lation parameters. However, neutrino oscillation results depend on the neutrino energy — a quantity 
which can not be measured directly but has to be reconstructed from the hadronic debris coming 
out of the neutrino-nucleus reaction inside the detector. A reliable reconstruction of the neutrino 
kinematics and the initial scattering process has to account for in-medium modifications and, in 
particular, for final state interactions inside the target nucleus. They can, e.g., through intra-nuclear 
rescattering, change particle multiplicities and also redistribute their energy. 

Those effects can be simulated with our fully coupled channel GiBUU transport model where 
the neutrino first interacts with a bound nucleon producing secondary particles which are then trans- 
ported out of the nucleus. We use a formalism that incorporates recent form factor parametriza- 
tions and apply, besides Fermi motion and Pauli blocking, important ingredients of the many-body 
problem such as mean-field potentials, in-medium spectral functions and RPA correlations. The 
modeling of final state interactions includes a large variety of possible interactions channels and, 
furthermore, particles with an in-medium width are transported off-shell. 

This article is structured in the following way: First we introduce our model for the interaction 
of neutrinos and electrons with bound nucleons. We then outline our transport model used to de- 
scribe final state interactions (FSI). Thereafter, some sample results for vA scattering are presented. 
Finally, we apply our model to describe recent observations of the LBL experiments MiniBooNE 
and K2K. 



2. GiBUU model 

Lepton induced scattering in the GiBUU model is treated as a two step process: First, the 
leptons scatter of nucleons embedded in the nuclear medium. Then, the outcome of this initial 
reaction is propagated through the nucleus, using a hadronic transport approach. More details can 
be found in Ref. 

2.1 Initial vertex 

We focus on the charged current (CC) (vA^ i~X) reaction, but discuss also the electromag- 
netic (EM) {£~N i^X) one used as a benchmark for our neutrino calculations. 

We treat the nucleus as a local Fermi gas of nucleons bound in a mean field potential. The 
total reaction rate for the scattering of a lepton with four-momentum k = {kQ,k) off a nucleon with 
momentum p = {E,p), going into a lepton with momentum k' = (^o,k') is given by an incoherent 
sum over all nucleons (impulse approximation) 

U'-'EM.CC _ / "*^EM,CC \ (2 1) 



d(0 dD. jr'i V /j 



with (0 = ko — k^, = —{k — k'Y and Q. = Z(k,k'). The cross sections on the rhs of Eq. (2.1) are 
medium-modified (see below). 
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Elementary input cross sections. In the intermediate energy region (^o ~ 0.5 — 2 GeV), the cross 
section is dominated by quasielastic (QE) scattering {eN e'N' and vA'^ £^N') and resonance 
excitation' (eN — > e'R and vA'^ i^R). Furthermore, we account for non-resonant single-pion 
backgrounds for both eN —>■ e'nN' and vA'^ — > i^nN'. Thus we assume 



dw dQ. dco dQ 



da^ dcj^'^ 
Vdwdn dcodQ.' 



(2.2) 



where da^^i^f, also contains contributions from resonance-background interference. 

Omitting phase space factors, the cross section for QE scattering and resonance excitation is 
given by [0] 



(2.3) 



d(0 dO. 

where p' = {E',p') is the four-momentum of the outgoing nucleon and £/ {E',p') gives the spectral 
function for the outgoing baryon. L^v is the leptonic tensor. 

The QE hadronic tensor Hq^ can be parametrized in terms of vector and axial form factors 
(see, e.g., our earlier work [||]). The vector form factors are taken from the latest analysis by Bodek 
et al. 10]; a dipole ansatz withM^ = 0.999 GeV [||] is used for the axial ones. 

The resonance hadronic tensor H^^ depends on the specific resonance. For spin 1/2 resonances 
with positive parity (e.g. Pii(1440)) we find for the hadronic current ^f/2+ 

(eV + 



^1/2 



</2-</2With 

Fp 



(2.4) 



and for states with negative parity (e.g. Sii(1535)) we use 7f/2_ 



For spin 3/2 resonances with positive parity as the P33(1232), we have j"'^ 



with 



/2J 

■XIX 

3/2+ 



3/2 ^3/2 



V. 



3/2 



+ 



Mr 



q-p 



q^p^)+g^^Cl 



(2.5) 



and 



"^3/2 



-q-f) 



+ 



(pA 



-q-p'^)+Ch''' + ^q-q^ 



(2.6) 



for the ones with negative parity (e.g. Di3(1535)) we use j"^ 



/3^2- — ^3/2 ~^3/2- approximation, 
resonances with spin greater than 3/2 are treated within the spin 3/2 formalism. 

The vector form factors F^2 (C^4 5(,^ present in CC scattering are related to the electro- 
magnetic transition form factors F^2 (^3^4 55) ^^^^ N = p,n and those again to helicity ampli- 
tudes ^, ||, ^ |l^, which can be extracted from electron scattering experiments. We apply the 



results of the recent MAID2005 analysis Ql IQ which includes 13 resonances with W <2 GeV; all 
of them are implemented in our model — a list is given in Table [H 

The lack of precise data renders the determination of the axial form factors difficult. Pion pole 
dominance and the PC AC hypothesis allow on one side to relate Fp to Fa (Cg to Cj ) and on the 



'mainly P33(1232) 
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name 


Mr [GeVl 


J 


I 


P 


[GeVl 

(J ' J 


p;riV/ptot 
/ 


axial coupling 




1.232 


3/2 


3/2 


+ 


0.118 


1.00 


1.17 


p., (1440) 


\A(yL 


1/2 


1/2 


+ 


0.391 


0.69 


-0.52 


Dn(1520) 


1.524 


3/2 


1/2 




0.124 


0.59 


-2.15 


Sii(1535) 


1.534 


1/2 


1/2 




0.151 


0.51 


-0.23 


S3i(1620) 


1.672 


1/2 


3/2 




0.154 


0.09 


-0.05 


Sii(1650) 


1.659 


1/2 


1/2 




0.173 


0.89 


-0.25 




1.676 


5/2 


1/2 




0.159 


0.47 


-1.38 


Fi5(1680) 


1.684 


5/2 


1/2 


+ 


0.139 


0.70 


0.43 


D33(1700) 


1.762 


3/2 


3/2 




0.599 


0.14 


-0.84 


Pi3(1720) 


1.717 


3/2 


1/2 


+ 


0.383 


0.13 


0.29 


F35(1905) 


1.881 


5/2 


3/2 


+ 


0.327 


0.12 


0.15 


P3i(1910) 


1.882 


1/2 


3/2 


+ 


0.239 


0.23 


-0.08 


F37(1950) 


1.945 


7/2 


3/2 


+ 


0.300 


0.38 


0.24 



Table 1: Properties of the resonances included in our model. The pole mass M^, spin J, isospin /, parity P, 
the vacuum total decay width Fq"^, the branching ratio into nN and the axial coupling are listed (see text for 
details on the extraction of the axial coupling). The resonance parameters are taken from Ref. 



other side to extract the axial coupling FpSS>) (^5(0)) (given in Table |l|). We assume a dipole 
form with M\ = 1 GeV for Fa and all except for the P33(1232). and are set to zero for all 
resonances except the P33(1232). 



For the A resonance, some experimental information is available from ANL [ |13| , |14( ] and 
BNL Applying the Adler model [|r|] where C^{Q^) = -C^{Q^)/4 and Cf (2^) = we can 
extract the dependence of Cj from these data. Any update of the vector form factors requires the 
axial ones to be refitted. Improving on the electromagnetic vector form factors without readjusting 
the axial ones [|8|] will result in a worse description of the neutrino data. Assuming that PCAC holds 
(i.e. the value of C5 (0) is unchanged^) and neglecting the small non-resonant background in the 
channel vp ^ jx^n^p, we find 



Cj{Q')=cUO) 



1 + 



b + Q^ 



1 + 



(2.7) 



with a 



-0.25 and b 



0.04 GeV^ and 



0.95 GeV (set "NEW"). Formerly, a = -1.21 and 



2GeV2 and 



1.28 GeV (set "OLD") were used [jlj, [14j, [IJ]. With the new parameters, 
good agreement with the ANL data is reached (solid line in Fig. |l|). In order to illustrate the 
sensitivity of the cross section to the refitted axial form factor, we show in Fig. |l| also the dashed 
curve which gives the cross section obtained with the new vector form factors but the old axial 
ones. 



^This is in contrast to the work of Hernandez et al. JT^ ] wiiich takes it as a free parameter. Since we rely on PCAC 
for all other resonance excitations (where no data are available) we prefer to keep it here, too. Furthermore, this coupling 
was extracted from the BNL data in Ref. |fl3] and found to be consistent with the PCAC prediction. 
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Figure 1: Differential cross section averaged over the ANL flux for two different sets of parameters 
describing Cj compared to the ANL data [jl^. An invariant mass cut W < 1 .4GeV is appHed. 



In the case of electro-production, there is a wealth of data which allows to determine the single 
pion background Ao^l^/AcoAQ.. This is done by subtracting the dominant resonance contribution 
from the total single-pion cross section [Q, |T^ 



dfi) dQ. dft) dQ. ^ do) dD. 



The total single-pion production cross section on the nucleon da^J /dcodQ. is taken from MAID [|TT]] 
Such a treatment is not possible in the neutrino case, since there are not enough experimental data 
to fix the additionally necessary six axial amplitudes, hence some simplifications are required^ and 
we assume 

dar? . da, 



where da^ is constraint by electron scattering data. The factor b^^ depends on the channel, vn — > 
l^pn^ (vp — > jX^K^p is assumed to be "background-free"): with b^^ = 3 and 
1.5 a reasonable agreement with the ANL data is reached as can be seen from Fig. ^ 



The full cross section Eq. ( |2.2| ) is shown in Fig. g for electron scattering off protons. The right 
peak is dominated by the A, the second and third resonance region are clearly visible (middle and 
left peak), the QE peak is not shown. The different contributions to the cross section are shown 
and compared to data. The inclusion of the Itt non-resonant background is necessary, but not 
sufficient to achieve a good description of the data at higher bombarding energies (right panel); 
multi-pion backgrounds should also be considered in the future. We further compare to the model 



of Rein and Sehgal [21 1, a model widely used in neutrino event generators, and find that this model 



underestimates significantly the electron data as observed also by Graczyk et al. [22]. In Fig. || 
we plot the neutrino induced total 71+ (left panel) and tt" (right panel) production cross sections; 
the different contributions are also indicated (full: resonances -i- background). The excitation of 
higher resonances is almost invisible in the isospin 3/2 channel (upper curves in the left panel, 
solid vs. long-dashed line). Again, we compare to Rein and Sehgal which gives for these integrated 
cross sections results similar to ours (comparing the dotted line to our resonance contribution). 



■'Besides the phenomenological ansatz we are following in this work, one can apply elementary models to estimate 
the single-pion non-resonant terms l^[l^, po| |. 
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Figure 2: Double differential cross sections for scattering of electrons off protons through resonance ex- 
citation and non-resonant processes as a function of the energy of the outgoing electron (sohd hnes). The 
dashed lines show the resonance contribution, and the dotted lines the outcome of the model of Rein and 
Sehgal. We compare to data from JLAB [H]. 
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Figure 3: Total (left) and tt" (right) production cross sections through resonance excitation and non- 
resonant processes as functions of the neutrino energy (solid lines) compared to the pion production data of 
of ANL (Refs. [|l|] (•) and [0 (■)) and BNL (x)). The different contributions are indicated by only 
resonances (short-dashed) and only A (long-dashed). The result obtained in the model of Rein and Sehgal is 
also shown (dotted). 



Medium modifications. The target nucleus is treated within a local Thomas-Fermi approxima- 
tion as a Fermi gas of nucleons bound by a mean-field potential Us^^ (P) r) which is parametrized as 
a sum of a Skyrme term depending only on density and a momentum-dependent contribution. 

The spectral function of a particle with four-momentum p = and mass M = ■\fp^ is 

given by 

'^^'^'P^ " ^(M2-M2-ReI(£,|)))2'+(ImE(£,p))2' ^^"^^^ 

with the self energy r(£',p) and the vacuum pole-mass Mq. It includes the effect of the momentum- 
dependent potential on the outgoing baryons and also accounts for the in-medium collisional broad- 
ening of the outgoing final states. We neglect the spectral functions of the initial states because their 
widths are considerably smaller than those of the outgoing nucleons [|2^. The imaginary part of 
the self energy is related to the full width, Ftot, in the medium, via ImE(£',p) = — MFtot, which is 
given by Ftot = TpB + rcoii- Due to Pauli blocking (PB) of the final state particles in the medium, 
the free decay width is lowered. On the other side, both the nucleons and the A resonances undergo 
collisions with the nucleons in the Fermi sea. This leads to a collisional broadening of the particle 
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width. To estimate this collisional broadening, we apply the low-density approximation 



rcoii(£,p)=E / a(£,p,p') v,eiPpB T^T^s, (2.11) 



where we integrate over all nucleon momenta in the Fermi sphere (FS). a(£', p, p') denotes the total 
cross section for the scattering of the outgoing nucleon/resonance with a nucleon of momentum 
p' in the vacuum; Vrei denotes the relative velocity of the particle and the nucleon, /^b is the 
Pauli blocking factor for the final state particles. The total cross sections are chosen according 
to the GiBUU collision term (for details see [|l|]). The real part of the self-energy is given as 
once-subtracted dispersion relation where the subtraction point is fixed by the mean fields. This 
procedure guarantees the proper normalization of the spectral functions 

In the nucleus, the elementary cross sections discussed above for resonance excitation and 
quasielastic scattering are evaluated with full in-medium kinematics accounting for the momen- 
tum-dependent mean field. Furthermore, also the flux and phase-space factors are evaluated with 
in-medium four-vectors. As an approximation, we use in the medium the same form-factor para- 
metrizations as in vacuum. Pauli blocking is taken into account by multiplying each cross section 



on the rhs of Eq. (gjj) with the Pauli-blocking factor. In particular, the momentum dependence of 
the potential and the collisional broadening of the nucleon improve the correspondence with the 
data considerably [Q]. 

2.2 Final state interactions 

The final-state interactions (FSI) of the produced particles are implemented by means of the 
coupled-channel semi-classical Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model 
[|l|]. Originally developed to describe heavy-ion collisions [25], it has been extended to describe the 



interactions of pions, real and virtual photons and neutrinos with nuclei |l^,^, 27, 28, 29, 3^]. 

In this model, we describe the space-time evolution of a many-particle system under the in- 
fluence of a mean-field potential and a collision term by a BUU equation for each particle species. 
A collision term accounts for changes (gain and loss) in the phase space density due to elastic 
and inelastic collisions between the particles, and also due to particle decays into other hadrons 
whenever it is allowed by Pauli blocking. The most relevant states for neutrino-induced reactions 
at intermediate energies are the nucleon, the A resonance and the pion. For the NN cross section 



and its angular dependence we use a fit to data from Ref. [|31|]. For the pion cross sections we use 
a resonance model with the background fitted to data as shown in detail in Ref. [29]. The decay of 
resonances into a pion nucleon pair is Pauli blocked if the momentum of the nucleon is below the 
Fermi momentum. We allow not only for the decay of the resonances, but also for the rescattering 
in the nuclear medium through processes like RN —>■ NN, RN —>■ R'N and, for the A resonance we 



also consider ANN NNN based on [32] 



In between the collisions, all particles (also resonances) are propagated in their mean-field 
potential according to their BUU equation. We emphasize again, that, due to rescattering effects 
in the medium, the nucleon and the resonances acquire an additional complex self energy leading 
to modified spectral functions obtained in a consistent way from the GiBUU cross sections (see 
Eq. ( |2.1 1[ )). Thus, the nucleon and the resonances are transported off-shell in our model. Thereby 
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Figure 4: Single proton knockout (left) and single tt" production (right) cross sections for v scattering off 
^^C (solid lines) through QE, resonance excitation and non-resonant processes. Possible origins, i.e. the 
initial processes, are indicated (QE by long-dashed, A by short-dashed lines). Note the different scale. 



we ensure that the particles are transported back to their vacuum spectral function when leaving 
the nucleus. 

In conclusion, FSI lead to absorption, charge exchange and redistribution of energy and mo- 
mentum as well as to the production of new particles. In our coupled-channel treatment of the FSI 
— in which the BUU equations are coupled through the collision term and, with less strength, also 
through the potentials — our model differs from standard Glauber approaches that do not allow for 
side-feeding and rescattering. 

Within the GiBUU model we performed for the first time a systematic study of how quasielas- 
tic scattering and resonance excitation are interconnected by FSI, focusing, on one side, on reso- 
nance induced nucleon knockout and "fake" CCQE events and, on the other side, on side-feeding 
effects in pion production [|3[ An example of the coupled channel effect in neutrino nucleus 
(here: ^^C) reactions is given in Fig. Q for p and Tl" yields which are of interest in current LBL ex- 
periments. We find, that even though the proton knockout is dominated by protons coming from an 
initial QE reaction, the secondary protons from an initial A excitation contribute significantly. FSI 
have also a considerable influence on CC neutrino induced pion production, where, due to isospin 
relations, 71+ dominate in the beginning. In the nucleus, they rescatter, get absorbed and undergo 
charge exchange reactions leading to a disproportionate population of the tt" channel. This has 



been earlier observed by Paschos et al. [33]. We stress that a correct understanding of the tt" yield 
is required for a correct identification of the neutrino flavor in LBL experiments. 

To conclude, we emphasize, that GiBUU is based on well-founded theoretical ingredients and 
it is the only model tested in various and very different reactions using the same physics input. 
In particular, an important prerequisite for any model for the interaction of neutrinos with nu- 
clei is, that it provides a good description of electron- or photon-induced reactions. Within the 
GiBUU model, extensive tests against existing data are possible and have been successfully per- 
formed [^, [l^, 34]. The GiBUU model is capable of incorporating the complexity of the nuclear 
many-body problem in an extensive open-source computer code which can be downloaded from 
our website [[I]]. 



3. Applications: MiniBooNE and K2K 

We shall now present some examples for the application of our model. 
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3.1 CCQE 

Charged current quasielastic (CCQE) events are commonly used in LBL experiments to de- 
termine the kinematics. Under the assumption that the nucleon is at rest within the nucleus, the 
neutrino energy has been reconstructed from QE events at the MiniBooNE experiment [|^] using 

^ _ 2{Mn - Eb)E^ - {El - IMnEb + ml) 
2[{MN-EB)-E^ + \k'\cosd^,] ' 

with a binding energy correction of Eb = 34MeV and the measured muon properties. With that, 
we obtain the reconstructed via 

Q^ = -ml + 2Ey{E^-\k'\cose^), (3.2) 

Two immediate questions are raised by this procedure: (1) How good is the identification of CCQE 
events? (2) How exact is the crucial assumption of two body kinematics for nucleons bound in a 
nucleus where many in-medium modifications are present? 

The experimental task is now to identify true CCQE events in the detector, i.e., muons origi- 
nating from an initial QE process. FSI might lead to misidentified events, e.g., an initial A whose 
decay pions are absorbed or which undergo "pion-less decay" can count then as CCQE event (we 
call this type of background events "fake CCQE" events). We denote every event which looks like 
a CCQE event by "CCQE-like". At MiniBooNE these are all the events where no pion is detected 
while at K2K these are all events where a single proton track is visible and at the same time no 
pions are detected. The two methods are compared in Fig. ^. The "true" CCQE events are denoted 
with the solid line, the CCQE-like events by the short-dashed one. Placing a cut only on pions, 
as MiniBooNE does, leads to a considerable amount of "fake" CCQE events (left panel, the short- 
dashed line is higher than the solid line). They are caused mainly by initial As via the mechanism 
described above; their contribution to the cross section is given by the dash-dotted line. On the 
contrary, less CCQE-like than true events are detected with the K2K method when one cuts both 
on pions and protons (left panel, difference between short-dashed and solid line). The final state 
interactions of the initial proton lead to secondary protons, or, via charge exchange to neutrons 
which are then not detected as CCQE-like any more {single proton track). We find that at K2K the 
amount of fake events in the CCQE-like sample is less than at MiniBooNE (compare difference 
between short-dashed and long-dashed line). We conclude that even if the additional cut on the 
proton helps to restrict the background, an error of about 25 % remains, since the measured CCQE 
cross section underestimates the true one by that amount. 

The flux averaged CCQE-like distribution for MiniBooNE is shown in Fig. |^ The influence 
of the fake CCQE events on the energy reconstruction is can be inferred from the left panel. The 
distribution obtained by reconstructing for the CCQE-like events via the formulas above is 
compared to the distribution of the true events (solid vs. dash-dotted line). At lower it is higher 
than the latter, but then it falls off faster. The difference between the two curves is caused by the 
fake CCQE events whose different muon kinematic affect the reconstruction. We also find that the 
reconstruction with the simplified formulas above turns out to be almost perfect when only true 
CCQE events (and not the whole CCQE-like sample) are taken into account. 

In the right panel of Fig. ^ we compare our calculation to the recent MiniBooNE findings: 
The MiniBooNE collaboration observed a dip at low compared to their standard Monte Carlo 
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Figure 5: Total QE cross section on '^C (solid line) compared to different methods on how to identify 
CCQE-like events in experiments (short-dashed lines). The contributions to the CCQE-like events are also 
classified (CCQE-like from initial QE (long-dashed), from initial A (dash-dotted)). 




[GeV^] [GeV^] 



Figure 6: Flux averaged distribution of CCQE-like events at MiniBooNE. The left panel shows the 
composition of the cross section (initial QE, initial A) and the reconstructed one. The right panel shows the 
influence of different in-medium modifications compared to the MiniBooNE tuned Monte Carlo output. For 
simplicity, only QE and A excitation were considered for the initial vertex. 



prediction [pSp. To reach agreement with the data, apart from changing the axial mass Ma by 
about 25 % with respect to values obtained in earlier experiments they had to modify Pauli 
blocking in their Monte Carlo description! Since MiniBooNE has not provided absolute cross 
sections, a direct comparison to data is not yet possible. Instead, we compare to their CCQE- 
hke Monte Carlo points^ which includes the aforementioned tuning of parameters and treat this as 
experimental result. We compare these points to our calculations including different in-medium 
modifications on the cross section. While the momentum dependent potential and the spectral 
function had significant influence on the double differential distributions [^, they are negligible 
here where we have integrated out one quantity, and do not improve the correspondence with the 
data. Polarization effects due to the strong interaction among nucleons modify the QE hadronic 



tensor. These are taken into account by including RPA correlations taken from Nieves et al. [ |36| ] 
which lowers the spectrum at the peak and leads to a good description of the MiniBooNE points 
without tuning the axial mass or any other parameter. We shall explore this in more detail in a 
forthcoming publication. 



by normalizing the area 
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Figure 7: Ratio of the total 7i;+ to the CCQE-like yield in K2K compared to data taken from Ref. [38|. 



3.2 CClTT 

LBL experiments demand for a realistic description of pion yields, for two main reasons: 
First, tt'^s are an important background in Vg appearance experiments which, through their photon 
decay, might affect the electron detection yielding also "fake" appearance events. Second, absorbed 
pions contribute to the CCQE-like background as discussed above. Therefore, a good description 
of neutrino induced pion production in nuclei taking into account complex FSI is necessary and 



oversimplifications are not justified. We have shown [37] that, in particular, the ANP model applied 
by Paschos et al. [|^] does not incorporate the well-known properties of the Nl^n dynamics in 
nuclei and, therefore, is not able to give reliable results for pion spectra in the energy region of 
interest. 

As an example, in Fig. ^ we compare the recent K2K data for the ratio of ;r+ to QE yields [38] 
with the output of the GiBUU model and obtain good agreement. 

3.3 Radiative A decay 

MiniBooNE finds in its — > Vg oscillation result an excess of electron like events for neutrino 
energies less than 475 MeV which is not yet understood [|^]. A possible source is the excitation 
of a A resonance via neutral current interaction followed by the radiative decay A yN. Since the 
MiniBooNE detector cannot distinguish between photons and electrons, this reactions gives rise 
to additional events in the low energy region. The major -induced background, however, are tt" 
coming from NC interactions detectable also via their photon decay products. Of particular interest 
for experiments is thus how the photon to TlP yield changes in the nuclear medium, depending on 
the A momentum and also the nuclear density. 

In the vacuum, a rough estimate gives 



-(A+/<') 



0.0056 



0.008, 



(3.3) 



where 0.0056 is the PDG branching fraction and 2/3 comes from the appropriate Clebsch-Gordon 
coefficient for A+Z" — > +N. In the medium — we took '^C as used in MiniBooNE — our calcu- 
lation has been performed as follows. First, A" and A+ resonances are set inside the nucleus with 
momentum and radius chosen randomly within a given range. Then, they are propagated out taking 
into account all kind of decays and collisions. Afterwards, we calculate the total and the photon 
cross section as function of the initial momentum (radius) of the A. We take into account only 
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Figure 8: Photon/ ^r" rate changes in medium as a function of A position (left) and momentum (right). 



those tt" which actually made it out of the nucleus after the final state interactions. With that, we 
obtain for the above ratio 0.019, which represents an increase of about a factor of 2.4 compared to 
the free case. 

Fig. |8] shows how the photon/Tr" rate changes in medium as a function of A momentum and 
position (solid: initial A*', dashed: initial A+). In addition, the vacuum estimate is shown by the 
long dashed line. In the momentum dependence one observes typical final state interaction effects: 
Slow As produce slow pions which are more likely to be absorbed in the medium than higher 
energetic ones which might pass through undisturbed. As expected, the medium modification is 
largest for those As which have been put in the middle of the carbon nucleus. However, one might 
expect, that the solid/dashed lines approach the vacuum value at a radius larger than the carbon 
radius. This does not happen here, because, as said before, we initialize the As at the beginning 
with a random momentum, therefore some of them can propagate into the nucleus and thus still 
undergo FSI which then again modify the spectrum. 

To conclude, the production of photons vs. tt" is enhanced in the nuclear medium due to 
complex pion final state interactions reflecting in a strong dependence of density and momentum. 



4. Summary 

We conclude that in-medium effects in vA scattering, and in particular FSI, are important for 
the interpretation of LBL oscillation experiments. The influence of nuclear many-body effects and 
final state interactions have to be treated with the same degree of sophistication as the primary 
production process. 
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